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1. INTRODUCTION 


Reflector antennas are widely used in communication satellite systems 
because of their relatively good radiation characteristics, low cost, and 
light weight. A central problem in the analysis of a reflector antenna is 
the secondary pattern computation. As sketched in Figure 1, for an inci- 
dent field from a feed located at , the problem is to calculate the scat- 
tered field E from a known reflector I at a far-field observation point r. 
Several methods exist for calculating the high-frequency asymptotic solu- 
tion of E (r), as explained below. 

(1) Physical Optics Method (PO, Fig. la) [1]— [7] . The induced 

A ^ . 

current on the reflector is approximated by 2n x H . An integra- 
tion of this current over the curved reflector Z gives the far 
field E S . 

(ii) Geometrical Theory of Diffraction (GTD, Fig. lb) [8] -[10]. At a 

> s 

far-field observation point r, the scattered field E consists of 

- - "f 

two terms: the reflected field on ray P^O , and the edge 

diffracted field on ray P^O^ . 

(iii) Aperture Integration Method (AI, Fig. lc) [1 ] , [ l 1 ]-[ 1 3] . The 
field on aperture plane Z is first calculated by tracing a 

cl 

reflected ray P^O P^ using geometrical optics theory and an edge 

diffracted ray P^O^P^ using GTD. Next, we integrate the field 

over Z via FFT to obtain the scattered far-field E S . 
a 

The accuracy of the above three methods is discussed below. As a 

■> jg 

reference, let us represent the exact solution of E by a high-frequency 
asymptotic series, namely, 
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£ s (r) ~ e [^ Q (r) + k 1 ^ 2 ^(r) + k ^(r) + •'•• ] * for k * ® (1.1) 

Then, we may summarize the accuracy and limitations of the three methods in 
the following table: 


TABLE 1 


ACCURACY AND LIMITATIONS OF PO, AI, AND GTD 


Methods 

Accuracy 

Limitation 

PO 

recover and partial A^ 

not accurate for 

AI 

wide-angle lobes 

GTD 

recover X and A, 
o 1 

predicts infinite 
field in main 



beam direction 
(caustics) 


In this report, we will study the main reflector far-field pattern 


using the aperture integration (AI) method for the following reasons: 


(i) Unlike the other two methods, AI gives the near-field (aperture 
field) as well as the far field. Most of today's large reflector 
measurements are done in a near field range. Thus, only AI 
provides a convenient theoretical check for the near field 
measurements. 

(ii) The accuracy of AI is comparable to the popular PO. The use of 
FFT in AI makes it numerically efficient. Furthermore, as will 
be discussed later, the present AI formulation is most suitable 
for extension to multiple (2 or more) reflectors. 

There exists an extensive list of published literature on AI - notably, 
Silver [1], Kauffman and Croswell HU, Acosta [12], and Hwang, Tsao, and 
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Han [13]. In comparison with these prior works, the present AI analysis may 
be considered as an extension in one or more of the following areas: 

(i) The surface of the reflector is completely arbitrary. It can be a 
numerically specified surface. 

(ii) The edge of the reflector is not restricted to a circular curve. 

It can be an arbitrary curve lying on an elliptical cone or 
cylinder. 

(iii) The divergence factor of the GO field is correctly computed. 

Hence, our analysis is not restricted to feeds located very close 
to the focal point (in which case the divergence factor is nearly 
unity and is ignored by several researchers). 

(iv) The edge diffracted field is included in the aperture field 

calculation. Near the incident and reflected shadow boundaries, 
two uniform theories £ 14 ] — [ 15 J are used so that the aperture 
field is continuous from the lit to the shadow region. 

In short, the present AI analysis represents a generalized and improved 
version of previous work. In particular, it is amendable for a convenient 
extension for analyzing multiple reflectors. 

The organization of this report is as follows: The description of the 

problem is described in Section 2. The' incident field, feed coordinates, 
and power radiated from the feed are covered in Section 3. Sections 4 and 
5 give a step-by-step procedure to compute the various contributions that 
make up the aperture field, namely, the geometrical optics and edge- 
diffracted fields. Section 6 covers the aperture field theory, the applica- 
tion of the FFT, and gain normalization. Numerical results and concluding 
remarks will be presented in Sections 7 and 8. 
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2. DESCRIPTION OF PROBLEM 


The geometry of the problem under consideration is sketched in 

V* 

Figures 2 and 3. A reflector £ is illuminated by the incident field from 
a point source at P^. The problem is (a) to determine the high-frequency 
asymptotic solution of the total field at an observation point on the aper- 
ture grid £ as shown in Figure 2 , and (b) to determine the secondary 
a . 

pattern using the FFT as depicted in Figure 3. In this section, we shall 
describe the various elements involved in the problem. 

2.1. Coordinate Systems and Time Convention 

The main coordinate system is the rectangular system (x,y,z), whose 
origin and orientation are arbitrarily chosen. In calculating the edge 
diffracted field which involves the boundary of the reflector, we employ 
a primed rectangular system (x"* ,y" ,z'' ), whose relation with (x,y,z) is 
explicitly stated later. The feed coordinate system (x^,y^,z^) required to 
describe the polarization and incident field is related to the main coor- 
dinate system by Eulerian angles [16}. This will be discussed in detail in 
Section 3. The time factor is exp(+jtut) and is suppressed throughout. 

2.2. Source 

We assume that the source has a well-defined “phase center" at P^ with 
coordinates (x-.y-.z^), and radiates a spherical wave denoted by (fi*,E*). 

If the feed is an array, it is necessary to consider each element in the 
array separately and superimpose their final scattered fields. 

2.3. Reflector £ 

The perfectly conducting surface is described by the equation 
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Figure 2. A reflector J with edge T being illuminated by 
the incident field from a point source at P^. 



Figure 3. Secondary pattern using the Fast Fourier Transform.. 


( 2 . 1 ) 


z = f(x,y) , for a < x < b and c < y < d 

It is not necessary to know the analytical form of the function f(x,y). In 

fact, the present computer program requires only a set of discrete data 

points ( x n ,y n ,f n ) with n = 1,2,...,N as the description of I. Those points 

are fitted by a cubic spline which gives automatically first and second 

2 2 2 2 2 

partial derivatives of f, namely, 3f/3x, 3f/3y, 3 f/3x , 3 f/3x3y, 3 f/3y . 
There are two requirements on the cubic-spline fit: (I) the data points 

can be distributed over a random grid, but they must be dense enough to 
describe the fine details of E; (ii) the domain of the data points 
(a < x < b,c < y < d) must be somewhat greater than the area defined by the 
boundary T of the reflector. Thus, we must know surface E in the shaded 
region in Figure 4 as well. Typically, the "width" of the shaded region is 
about 3 to 4 wavelengths. Our present program contains an extrapolation 
subroutine, which automatically extends E outward “smoothly to obtain the 
necessary data points in the shaded. region. " The final scattered field, 
for all practical purposes, is independent of the surface outside T. 
Consequently, the exact manner in which the extrapolation is done is 
unimportant. 

2.4. Boundary T 

Two types of reflector boundaries are most frequently used in prac- 
tice, and they receive our special attention. 

(a) Cylinder Case. In the first case, T is the intersection of sur- 
face E and an elliptical cylinder (Figure 5a). The parameters of the 
cylinder are: 
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Figure 5. Two examples of boundary T of the reflector. 
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(x^ ,y c ) = center of the ellipse 
(K 1 ,K 2 ) = semiaxis along (x,y) direction. 

Using a point on the axis of the cylinder as the origin, we introduce the 
second rectangular coordinates system (x',y',z') such that 

x' = X “ X 

c 

< y' = y - y c (2.2) 

z' = z 


The curve T may be described by a parametric equation with parameter <\>' , 
0 < < 2 ir, 


x' = gjU") = p' cos y 


r: 


g 2 (<j>') = p' sin 


{ z ' = g 3 ((f>') = f ( x = x c + g [( y = y c + g 2 ) 


(2.3) 


where 


P = 



COS <t>' 

2 

4 _ 

sin (j>" 

2 


K 1 





- 1/2 


(2.4) 


The projection of T on the x-y plane is always an ellipse. 

(b) Cone Case. In the second case, T is the intersection of surface 
Z and an elliptical cone (Figure 5b). The axis of the cone lies in the 
y' - z' plane and its parameters are 

(x = 0 ,y = 0 ,z = -p) = tip of cone 

0 ^ = inclination angle of cone axis measured from z-axis 
( 0 i , 83 ) = half-cone angles in the x' - z' and y" - z' plane 
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Using the cone tip as the origin, the primed coordinates (x - * ,y' ,z') are 
introduced such that 

x" = x 

y' = -y cos 0^ - (z - p) sin 0^ 

z' = y sin 0 3 - (z - p) cos 0 3 (2.5) 

Using parameters <j>', the curve T is described by 
x' “ = P' cos <f>' 

y' = g 2 ( ♦') = P' sin <p" 

z" = g 3 ( <*>') = p'[(cot 0 cos <r ) 2 + (cot 0 sin <f>') 2 ] l/2 (2.6) 

To determine p' as a function of we must solve the following nonlinear 
equation: 

. f(x,y) - p + y' sin © 3 + z' cos © 3 = 0 (2.7) 

For a given <j>^, there is a unique root p' from (2.7). The pair (^-.p') 
gives the desired relation p' = p'($') t which is fitted by spline functions. 
The projection of F on the x-y plane is, in general, a pear-shaped curve. 

(c) Arbitrary Case. In addition to the above two frequently used 
special cases, F may be an arbitrary curve described by 

x' = gjUO.y' = g 2 (<l>'),z' = g 3 (<K) 

where (x' ,y" ,z') is related to (x,y,z) by either (2.2) or (2.5). The func- 
tions (gj,g 2> g 3 ) can be specified either analytically or numerically. 
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2.5. Aperture Grid Points 

A general point on the aperture plane E is denoted by P 9 with coor- 
dinates (x^, y 2 ,z 2 ). Each point on E a as shown in Figure 3 is fed into the 
computer for repeated calculations of the scattered field. It is essential 
that successive observation points are adjacent, because of the following 
fact. In determining the reflection (specular) points on the reflector, we 
make an exhaustive search only for the first observation point in a batch. 
From the second point on, we use the reflection point of the previous 
observation point as the initial guess for the current reflection point. 

It is only when successive observation points are adjacent that such an 
initial guess ensures fast convergence. In the cases that were considered, 
only one iteration was needed to obtain the reflection point for all obser- 
vation points other than the first observation point. 

2.6. Method of Solution 

+ £ -+i 

For a given incident field (H ,E ) from the source at Pp the asymp- 
totic solution of the total field (ti^.E*") at point P 2 is determined by 
Keller's geometrical theory of diffraction (GTD) [17], Explicitly, the 
total magnetic field is asymptotically given by 

Keller: H t (P 2 > = H g + H d + 0(k _1 ), k > » (2.8) 

q 

Here H , the geometrical optics field, is of order k relative to the inci- 
dent field and is the dominant term. The second term H d is the edge- 
diffracted field and is of order k ^ 2 . It is well-known that H d becomes 
infinite and (2.8) fails if observation point P 2 is close to the incident 
or reflected shadow boundary. In the latter case, we will use the uniform 
asymptotic theory (UAT) [14], (18]— [21] , which amounts to replacing H d in 
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(2.8) by a new terra H° (with capital D): 

UAT: H C (P 2 ) = H g + H° + 0(k _1 ) , k > - (2.9) 

Once H is found, we calculate E from it by using the fact that H & , H , 
and are all the so-called "ray fields" which are locally plane waves. 

Once the tangential fields (E or H ) on I are obtained, the Fourier 

cl 

transform of these fields will essentially produce the secondary pattern. 
This will be discussed in detail in Section 6. 
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3. INCIDENT FIELD 


The definition of the incident field is given. The incident power 
from an arbitrarily polarized feed is derived so that the secondary pattern 
has the correct gain level. Finally, the computation of the incident field 
at a point on the reflector in the reflector system is described. 

3.1. Definition of Incident Field 

The surface current at the radiating aperture of the feed element may 
be expressed as 

l(x,y) = I(x ae^ + yb) (3.1) 

where (a,b,<|>) are real and 

a 2 + b 2 =■ 1 . (3.2) 

By choosing (a,b,\|>), one may obtain any feed polarization. Table 2 shows 
the values of a, b, <j> corresponding to linear and circular polarizations. 


TABLE 2 

VARIOUS FEED POLARIZATIONS 

a b tj; 


linear x 

1 

0 

0 

linear y 

0 

l 

0 

RHCP 

\/7T 

1 / 7i 

90° 

LHCP 

i/Ji 

1 /7I 

-90° 
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The radiated electric field due to 5 given by (3.1) is 

E ~ f (0 ,<*,) (3.3) 

where f(0,<j>) is the active element pattern. The function f(0,<j>) may be 
approximately expressed by 


t(0,<}>) = 0U E (0)(ae^ cos $ + b sin <J> ) + $U H (0)(b cos <J> - ae^ sin <t> ) (3.4) 

where ■ 

Ug(0) = E-plane active pattern 
U H (0) = H-plane active pattern 

Typically, these active patterns may be approximated by (cos 0)^, i.e. , 

(3.5a) 
(3.5b) 


From (3.3) and (3.4), the radiated electric field is given by 


U E (0) = (cos 0) qE 
U H (0) = (cos 0) qH 

3.2. Incident Power Radiated 


^-jkr 


- — [U E (0)(ae^ cos $ + b sin <j> )9 + U H (0)(b cos $ - ae^ sin <}>)<)>] 




The power radiated, assuming forward radiation only, is 
2ir it / 2 * > 

p = f / r 2 sin 0d0d$ 

inc i-0 0=0 Z o 


(3.6) 


(3.7) 


where 


Z = 120rr ohms, 
o 


Using (3.2), (3.6), and (3.8), the incident power radiated is 


(3.8) 


tt/2 


P inc = tW / i U E U E + U H U h] sln 0d9 


(3.9) 
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(3.10) 


For (cos 9) q type patterns, 


q E + q H + 1 


inc 60(2q_ + l)(2q u + 1) 

h n 


3.3. Determining the Incident Field on the Reflector 

The geometry of this problem is illustrated in Figure 6. The feed 
coordinate system (x f ,y f ,z f ) and the main reflector system (x,y,z) are 
related by Eulerian angles Yj.YjtY-j* In the most general case, these 
Cartesian systems can be aligned by three rotations. The angles of these 
rotations are known as Eulerian angles. Figure 7 illustrates the Eulerian 
angles The definitions of these angles [22] are as follows. 

Angle Yj describes a counterclockwise (ccw) rotation about the z axis which 

brings the x axis to the x” axis aligned with the line of nodes (line of 

intersection between xy and x^y^. planes), angle defines a rotation about 

the line of nodes in a ccw sense as indicated so that this brings the z 

axis to z^, and angle Y^ is another rotation about the z^ axis and aligns 

A A 

the x" axis with the x f axis in a. ccw sense. Typically, x f and z f are 
given. Let these unit vectors be expressed by 


Xf = *x + x^y + x 3 z 


z f = z^x + z^y + z 3 z 


(3.11) 

(3.12) 


Then, the Eulerian angles are given by 


Yj = tan 


-1 


'2 J 


(3.13) 


Y o = tan 


-1 


/ 


1 - z„ 


(3.14) 
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y 3 = can 


-1 


( z l x 2 “ x l z 2 J 


(3.15) 


From [16], it can be shown that 


AAA 


r lT = r A * *1 T 

[Xf y f Z f ] = A[x y zj 


(3.16) 


where 


A = 


cosy ^ cosy ^-siny^ cosy 2 siny ^ siny ^ cosy 2 +cosy ^ cosy .^siny^ siny^irtYj 

-cosy^siny^-siny 1 cosy 2 cosy ^ cosy ^cosy 2 cosy 2 _ siny ^siny 2 siny 2 cosY 2 

-cosy ^siny 2 


siny ^ siny 2 


cosy , 


(3.17) 

Thus, a point with coordinates (x,y,z) in the main coordinate system on 

F F F 

E has coordinates (x ,y ,z ) in the feed coordinate system given by 


F 


' 

X 


X - Xj 

F 

= 


y 

= A 

y - yj 

F 



z 


Z - Zj 


(3.18) 


The corresponding spherical coordinates are 


r , f\ 2 j. t *\2 ^ , F. 2-i 1/2 
r = [ (x ) + (y ) + (z ) j 


0 f = cos 1 (z F /r f ) 


- ~ l r F / F, 
= tan (y /x ) 


(3.19) 
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From the feed function given by (3.6), the incident E-field may be 
obtained. To find the incident H-field, the following equation is used: 

* > 

5 i„cMH r \ (3.20) 

r f r o 

where i s given by (3.8). Next, the incident field is converted from 
spherical components to rectangular components using 


V 


sin8 £COS<j> £ 

COS0£COS4)£ 

-s in<£ £ 


n 

>■*» 

i 

H. 

y f 

3C 

sin0£sin<}>£ 

cos0 £sin}> £ 

COS$ £ 


\ 

H 

L z f J 


COS0 £ 

-s in8 ^ 

0 


-V 


(3.21) 


Finally, the incident H-field in the (x f> y f ,z f ) system is 
(x,y,z) system using the following equation 


'h 


'h 

X 

=T 

X f 

H 

= A 1 

H 

y 


y f 

H 


H 

_ Z 


L Z f J 


converted to the 


(3.22) 


where A is given by (3.17). 
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4. GEOMETRICAL OPTICS FIELD 


The geometrical optics field H® in (2.8) consists of two components: 
the incident field if*" and the reflected field i? r which is calculated in 
this section. It should be noted that the incident field at an observation 
point on the aperture grid is taken to be zero because the incident field 
does not contribute to the secondary pattern. 

4.1. Reflection Point 

For a given source point P^ and an observation point P 2 (Figure 2), ^ 
reflection point 0 r may exist on the reflector E, and we denote its coor- 
dinates by (x,y,z=f (x,y)). The vectors 

= x(x - Xj) + y(y - y^ + z(f(x,y) - (4.1a) 

? 2 = x( x 2 ~ x) + y(y 2 - y) + z(z 2 - f(x,y)) (4.1b) 

r r 

are the connecting vectors between P^ and 0 , and 0 and P 2> respectively. 
The condition on the reflection point is that the distance (d^ + d 2 ) must 
be stationary, i.e. , 

(dj + d 2 ) - 0 , (d t + d 2 ) - 0 (4.2) 

which is explicitly given by 

-J- {(x - Xj) + [f(x,y) - z x ] -||} + {(x - x 2 ) + [f (x,y) - z 2 ] = ° 

■ 

{(y - y x > + [ f (x,y) - z i 1 17} + {<y - y 2 ) + [ f < x »y> ~ z 2 i 17} = 0 

20 


(4.3a) 



A root (x,y,z=f) of the two nonlinear equations in (4.3a) gives the loca- 
tion of a reflection point. For a given and P 2 , there may be none, one, 
or more than one reflection point. It may be shown that (4.3a) is equiva- 
lent to the .satisfaction of Snell's law. 

r 

The system of equations (4.3a) can be also satisfied if , 0 , and P^ 
are collinear. Such a spurious root may be eliminated by an additional 
condition 


r 1 

x - X. x - X- 

. * 4 . ^ 

2 

+ 

y ~ y l , y “ y 2 ] 

2 

• *1 
z- - z; z - zy 

- + - - - 

d l d 2 J 

1 

l d l d 2 J 

f 

d l d 2 J 


(4.3b) 

where 6 is a small positive number. We set 6 = 0.0001. 

A root of (4.3) may or may not fall inside the boundary T of the 
reflector (Figure 4). Thus, for each root (x,y,z=f) or its corresponding 
coordinates ( x" , y" , z ' ) in the primed system, the following test must be 
performed. If 

/OO 2 + (y') 2 1 Jl 1 (4.4) 


then the root is inside T and it is indeed a reflection point on the 
reflector. If (4.4) is not satisfied, then the root should be discarded. 
The parameter 2, ^ in (4.4) is given by 


L = [(x ') 2 + (y'/j 


+ \2 


1/2 


•)2 


•12 


- 1/2 


w 


K„ 


if T lies on an elliptical cylinder (Figure. 5a); and 


(4.5a) 



2 ' 2 1/2 

|z' \[(x) + (y ') ] [(x' 


cot 0 j )‘ 


+ (y' cot 8-)“] 


- 1/2 


(4.5b) 


if T lies on an elliptical cone (Figure 5b). 
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4.2. Formula for Reflected Field 


The reflected magnetic field at is given by 

H r (P 2 ) = (DF) e jkd 2{H 1 (0 r ) - 2[ H i (o r )- N] N} (4.6) 

-►i 

which is given in terms of the incident field H at the reflection point 

r r 

0 , the surface normal N of the reflector at 0 , and a divergence factor 

A A a 

DF. We choose N pointing toward the source; thus, (N»z) is always greater 

A 

than zero (Figure 8). Explicitly, N is given by 


N=A(-fx-fy+z) 
x y 


(4.7) 


2 2 - 1/2 

where A = +(f^ + f + 1) and the subscript x of f^, for example, means 
partial derivative with respect to x. The divergence factor in (4.6) is 


DF 


A + (d 2 /Rj V A + (d 2 /R 2 ) 


(4.8) 


where the square roots take positive real, negative imaginary, or zero 
value (so that DF is positive real, positive imaginary or infinite). 
(Rj.Rj) are the principal radii, of curvature of the reflected wavefront 
passing through 0 C . Their computation is given next. 


4.3. Curvatures of Reflected Wavefront 

We use the formulas given in Section IV of [18] for calculating 
(Rj,R^). The three orthonormal base vectors of the incident pencil are 
chosen to be (Figure 8) 


i 



A A ^ A A 

y * * 3 x(z - Zj) - z(x - Xj) 

|y x *3! [( z - Zl ) 2 + (x - x : ) 2 ] 1/2 


(4.9a) 
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(4.9b) 


X 2 


X Xj 

1*3 * 


■'i 

X 3 = 


x(x - x^ + y(y - y^ ) + z(z - 2 ^) 

7 7 N 2 7 ,2 , , 2 , 1/2 

[(x - xp + (y - yj) + (z - Zj) ] 


(4.9c) 


r 

where (x,y,z) are the coordinates of the reflection point 0 . Those of the 
reflected pencil are chosen to be 


A A 1 A • A A 

x^ = Xj - 2(xj»N)N 


(4.10a) 


A A • A * A A 

x^ = x 2 - 2 (x 2 *N)N 


(4.10b) 


- r x(x 2 - x) + y(y 2 - y) + z(z 2 - z) 

x 3 = 77 ~2~— 2 7 7 s2, 1/2 

l(x 2 - x) + (y 2 - y) + (z 2 - z) ] 


(4.10c) 


A J * A J * A f 

Note that (4.10) chosen above is a left-handed system, i.e. , x j x x 2 = -x^. 

r it 

This choice, of course, does not affect the final solutions of (R^,R 2 ). 

The three orthonormal base vectors of reflector 2 at O r are chosen to be 
x + zf 

(1 + fV /2 (4.11a) 


''Z 

X 1 = 


y + zf 


x„ = 


(1 + f 2 ) 1/2 

y 


(4.11b) 


= N 


(4.11c) 


From (4.9) and (4.11) the elements 
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P 1 = x i *x J ' , m,n = 1,2,3 
mn m a ’ * * ’ 


(4.12) 


can be calculated with the results 


(z “ Zj) ~ f x (x - x L ) 


11 [1 + 4 ]1/2 [(x " x i )2 + ( z " > 2 1 1/2 


(4.13a) 


-f (x - X, ) 

y r_ 


12 M , , 2 , 1/2 .. .2 . . . 2 , 1/2 

[1 + f ] [ ( x - x L ) + (z - ] 


(4.13b) 
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-(x-x 1 )(y-y 1 )-f x (y-y 1 )(z-z 1 ) 

[ l+f x J {(x-x 1 ) 2 (y-y 1 ) 2 +[(z-z 1 ) 2 +(x-x 1 ) 2 ] 2 +(y-y 1 ) 2 (z-z 1 ) 2 ) 1 ^ 2 


(4.13c) 
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■ (z-z 1 ) 2 +(x-x 1 ) 2 -f y (y-y 1 )(z-z 1 ) 

[ 1+f 2 ] 1 ^ 2 {(x-x 1 ) 2 (y-y 1 ) 2 +[(z-z 1 ) 2 +(x-x 1 ) 2 ] 2 +(y-y 1 ) 2 (z-z 1 ) 2 } 1 2 


(4.13d) 


P 33 = 17 A[f x (x " x l> + f y (y ' y l } _ (z “ z l )] 


(4. 13e) 


The first four elements in (4.13a) through (4.13d) form the 2x2 matrix 

= ~ it = i 

P . Because of the particular choice in (4.10), we have P = P . The 

curvature matrix of the incident pencil is 


Q 1 - 


‘I 1 

0 


0 

h l 


(4.14) 
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r 

The curvature matrix of reflector E at 0 is 


Q E - 


A 2 (eG - fF) A 2 (fE - eF) 


|_A 2 (fG - gF) A 2 (gE - fF) J 


(4.15) 


where 


E = 1 + f 2 , F = f f , G = 1 + f 2 
x x y y 


e = -Af , f = -Af , g = Af 

xx xy ■ yy 


A is defined just below (4.7) 


= 1T 

The desired curvature matrix Q may be calculated from the following matrix 
equation 


Q r ='. Q 1 + ,2 P 3 3 [(P i ) T ] 1 Q E (P i ) 1 


(4.16) 


Let us denote the four elements of Q r by 


Q r = 


l Q U Q 12 


P21 Q 22 


(4.17) 


Then the desired radii of curvature of the reflected wavefront at 0 r are 
given by 


_r ’ „r 2 ^ Q ll + Q 22 ) ± ^ (Q 11 + Q 22 J - “ 4(Q 11 Q 22 Q 12 Q 21^ 


R 2 


(4.18) 
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r r r 

Both and are real. Their signs have the following meaning: If R^ 

is positive (negative), the corresponding normal section of the reflected 

wavefront is divergent (convergent). The same convention applies to R^. 

Two final remarks about the calculation of the reflected field: 

(i) For a given and ? 2 > there may be more than one reflection point on 
E. Then the total reflected field is the superposition of the contribu- 
tions from each reflection point. If there is no reflection point on £, 
the reflected field is zero; (ii) If the reflection point is close to the 
boundary T, we still calculate its reflected field in the usual manner. We 
shall adjust the diffracted field H later by using UAT so that the total 
field H C in this case is correct. 
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5. EDGE-DIFFRACTED FIELD 


5.1. Diffraction Points 

To calculate the edge-diffracted field, the diffraction points on the 
boundary T of the reflector must be located first. Consider a source point 
at (x'.y'.zp and an observation point ?2 at (x£,y£,z£) , with their 
coordinates given in the primed system (Figure 9). A diffraction point 
0 d with coordinates (x' ,y' ,z') can be determined from the law of diffrac- 
tion 

7 - I «t = -J- £ *t (5.1) 

d 3 3 d 4 4 

Here t is the unit tangent of T at 0 d , and 

d 3 = x ”"( x - xp + y'(y' - yp + z'(z' - zp (5.2a) 

d 4 = x-(x- - x') + y'(y' - y') + z'(z' - z') (5.2b) 

From Fermat's principle, (5.1) is equivalent to 

ir <d 3 + V - 0 < 5 - 3 > 

From (5.2), (5.3) may be written as a nonlinear equation for unknown $■*: 


g l - *1 . g l - X 2 


g 2 ~ y{ g 2 “ y 2 

ag 2 

d 3 d 4 

3 y 

d 3 d 4 

. 4 

Ty 



(5.4) 


where (gj,g 2 >g 3 ) are defined in (2.3) if T lies on a cylinder, and in (2.6) 
if F lies on a cone. A root of of (5.4) determines a diffraction point 
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on F. Depending on the geometry, there exist examples where as many as 
four diffraction points have been found. 


5.2. Formula for Diffracted Field 

Corresponding to each diffraction point 0 d , there is a contribution 
to the diffracted field H d in (2.8). Following Equation (5.21) of [18], the 
formula for such a contribution reads 


H d (P_) = g(kd,j 


1 


1 


[3DV + aDV] 


4 /I + (d 4 /Rj) sin 3 


(5.5) 


Here g is a cylindrical wave factor 


g(x) = 


2/ 2m 


exp 


- j x + T 


(5.6) 


The other factors used in (5.5) are explained below. 


5.3. Divergence Factor 

The square root in (5.5) as usual takes positive real, negative imagi- 
nary, or zero value. is a radius of curvature of the diffracted 

wavefront passing through 0 d . It may be calculated from Equation (5.11) of 
[18] which reads 


L . + -JE_ 

d 3 sin 2 0 


!_5 J_ 

d 3 3 d 4 


^4 * n 


Here 8 is the angle between tangent t and d^ (Figure 9), 
is given by 


sin 8 


= \/l - Q 2 


(5.7) 

The factor sin 8 
(5.8a) 
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where 


o - pH<*' - - y i ) ir + <z '- z '0w 

Yag ] 2 fa g 1 2 fag ] 2] 1/2 
3^j a*- a*' 


(5.8b) 


(5.8c) 


All the derivatives In (5.8) and in the remainder of this section are eval- 
uated at the diffraction point 0^, whose coordinates are The 

curvature ic and normal n of the curve T at 0 ^ can be calculated from 


Equation (13.9) of [23]. The results are 


(5.9a) 


n PR g 2 g 3 


> > • 


g 3 g i ■ g i g r 


§2^ j 


(5.9b) 


where the prime on g' signifies the partial derivative with respect to <b' 

n 


8 = t(z'g' 3 ' - 83 K 2') 2 + U' 3 s\' - K[fT 3 ') 2 + )2] 1/2 (5.10) 

Summarizing the results in (5.7) through (5.10), we obtain the final 
expression for the divergence factor of the diffracted field 


DFD = 


/I + (d./R,) /I + G 
4 1 


(5.11) 
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where 


G = __ + 


^3 P 4 (sin e) 2 


x' - x' x' - x' 

L + _ 


lg' 3 (z' 3 g''-g'g' 3 ') - g' 2 (g' 1 g' 2 '-g' 2 g'') 


r-r,r- y- 

d 3 \ 


) 


lg\(g\g' 2 -g 2 g\') ~ g' 3 (g 2 g' 3 -g 3 g 2 ')] 


z' - z' z' - z' 
1 + z 


d 3 d 4 


lg' 2 (g' 2 g 3 ' - g' 3 g 2 ') ~ g\(g' 3 g\' ~ g\g' 3 ')]j 


5.4. Diffraction Coefficients 

s H 

The soft and hard Keller's diffraction coefficients D and D are 
defined in Equation (5.22) of [18], namely, 

D s > h - X 1 * X r 

= -sec ($ - $*) ± sec i. + $*) (5.12) 

The angles <J> 1 and <j> are shown in Figure 10. Because of the fact that 

A A A A 

N*z > 0, it can be shown that vector t x N is tangent to Z at 0 r , and 
points away from (not toward) E. We calculate <j>* and from the relations 
(Proj <L). (t x N) 

cos 4. 1 = (-1) (5.13) 

I Proj ^ 3 | 

(Proj d 4 )*(t x N) 
cos <{> = 

I Pr °j d 4 | 

Here Proj d^ is the projection of d^ on the plane perpendicular to 
namely, 


(5.14) 


t. 
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(5.15) 


4- -V + 

Proj d 3 = d 3 - (d 3 »t)t 


x S 31 + y"S 32 + z S 33 


where. 


S 31 


d 3 Q 3g l 
P 


d 3 Q 3g 2 

^32 ^ y “ ~ a*' 


1' P 3<t>' 


d 3 Q 3g 3 


S 33 (z> Z V P 3<p r 


Similarly Proj 5^ is given by 

Proj d^ = x + y S^ 2 + z S^ 3 (5.16) 

where 


d 4 Q a«i 

5 41 = ( x 2 ~ x ' ) p~ 'W 

d 4 Q 3g 2 

5 42 = (y 2 " y " ) " “777 

d<J> 

d 4 Q 3g 3 

5 43 = (z 2 - z '> “ 

A j 

The normal N of the reflector at 0 U is given in (4.7). For the present 
application, it is convenient to change its base vectors to those in the 
primed coordinate system. The result is 


N = x'Nj + y'N 2 + z'N 3 (5.17) 


where 
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A 


(5.18) 


N 1 " _Af x • N 2 " - M y * N 3 


If boundary T lies on an elliptical cylinder; and 


N 1 = 


N 2 = -A(-f y cos 0 3 + sin 0 3 ) 


N_ = -A(f sin 0 - cos 0.) 
3 y 3 3 


if T lies on an elliptical cone. Substitution of (5.15) through (' 
into (5.13) and (5.14) gives the final formulas for calculating 

cos - [SjjH, + S 32 M 2 + S 33 M 3 nS 31 + S 32 + S‘ 3 J 


2 2 2 - 1/2 

cos = [s 41 Mj + s 42 m 2 + s 43 m 3 ](s 41 + s 42 + s 43 ] 


where 


M =^- 


N, 


Jh 

2 dp' 


- N. 


3g 2 

3 Ihp 


M„ = 


3g, 3g 3 

L n 3 W " N i TV. 


M 3 " P 


3g 9 ag. 

m — £ _ m — L 
N 1 dp' dp' 


The solutions of <j> and <J> are subject to the following tests: 

(i) (J> i takes the value in the range (0,tt) if T 3 >_ 0, and the 
range ( 7r , 2 vr) if T 3 < 0, where 
T 3 = (-Proj ct 3 ) -N = “S 31 N 1 - S 32 N 2 - S 33 N 3 


(5.19a) 

(5.19b) 

(5.19c) 

.19) 
and <j>: 

(5.20) 

(5.21) 


(5.22) 
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(ii) <j> takes the value in the range (0,ir ) if VO, and the range 
( 7T , 2u ) if < 0, where 


T 4 = (Proj d 4 )»N - S 41 N l + S 42 N 2 + S 43 N 3 


(5.23) 


(iii) If the observation point P 2 is exactly on the incident shadow 
boundary, then 


<j> - <j) - tt = 0 


(5.24) 


(iv) If the observation point P 2 is exactly on the reflected shadow 
boundary, then 


$ + <}> - u = 0 


(5.25) 


When the observation point is exactly on the incident shadow boundary, the 
first factor in (5.12) becomes infinite, and causes computational dif- 
ficulty. A simple remedy is to shift P 2 slightly whenever (5.24) is 
satisfied. It should be remarked that when P 2 is near but not exactly on 
the incident shadow boundary, x 1 is large but finite. It does not cause 
any computational difficulty at the moment. Later on, the diffracted field 
in this case will be modified to become by using the uniform asymp- 

. "► £ r 

totic theory so that the total field H is correct. Similarly, x in 
(5.12) becomes infinite if the observation point ? 2 falls exactly on the 
reflected shadow boundary. Hence, we shift P 2 slightly when (5.25) is 
satisfied. 


5.5. Spherical Components of Incident Field 
+ i +i 

Fields Hg and in (5.5) are the two spherical components in the 

*i ~± 

directions of Q and a of the Incident magnetic field H evaluated at 
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The 


Then 


Here 

They 

if r 
f 


if r 


*i 

base vectors (3 ,a ) are shown in Figure 11 and may be calculated from 


A 1 A -t + 

3 - (o x d 3 )/d 3 

a 1 = (t x d 3 )/(d 3 sin 3) 


(5.26a) 

(5.26b) 


it may be shown that 


h 3 - inrj t5i( ° d) *‘ ] 


-i 


P sin 3 




3 go 


ago 


H 1 3$" + H 2 3$' + H 3 


(5.27a) 


H 1 = 


a Pd 3 sin 3 



3 §2 

3go 

H 1 

_ar 

- Z P - 3*' 


+ H„ 


3go 3g. 

w (K ' ~W (z ' - ° 


+ H. 


3g. 3go 

a y ( y' ~ yp - (x ' " 


(5.27b) 


AAA 


(H^,H 3> H 3 ) are the components of H in the directions of (x',y',z'). 
are given by 


H 1 = H x ( ° d) » H 2 = H y (0<1) * H 3 = H z ( ° d) 


(5.28a) 


lies on elliptical cylinder; and 


H, = 


H*(0 d ) 


H„ = 


-cos 0. H ± (0 d ) - sin 0- H i (0 d ) 
3 y 3 z 


H_ = sin 0, H i (0 d ) - cos 0- H i (0 d ) 
J j y _> z 


(5.28b) 


lies on elliptical cone. 
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Figure 11. Spherical base vectors (fl*, a*) and (0, a). 
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5.6. Rectangular Components of Diffracted Field 


The diffracted field H calculated from (5.5) must be eventually 

AAA 

expressed in terms of rectangular components in the directions of (x,y,z), 
so that it can be conveniently superimposed with the geometrical optics 

^cr * A 

field H in Section A. In (5.5), the two spherical base vectors (8 ,a) 
(Figure 11) may be calculated from 


3 - (et x d 4 )/d 4 (5.29a) 

o = (t x 2 )/(d sin 8) (5.29b) 


The diffracted field expressed in terms of the primed base vectors are 


H d (P 2 ) = x"hj- + y'h 2 + z'h 3 


(5.30) 


where 


n 8 


a, = 


1 Pd. sin 8 
A 


= 


2 Pd. sin 8 
A 


= 


3 Pd. sin 8 
A 


8 1 d, [a 3 (y 
4 


s 2 ‘ 77 <V 2 ' - 


+ D S 

a 

n 

H 1 ] g(kd ) (DFD) 
a A 

i 

sin 8 



3g 3 


- 

T* 

i 

if 

- y' 2 > - 3^ 

(z' 

- z') 

3 *1 


3g 3 


1 

3<t>- 

( Z “ 

Z V 3$' 

(x' 

-x 2 ) 



3g l 


- 

3<ji' 

( x* - 

+ \ A 

~ x 2 ; “ 3 <J> " 

(y' 

- y' 2 ) 

-yp 

- 

a 2 (z' - z' 2 )] 



- z p 

- 

a 3 (x" - x 2 )] 




n - 1,2,3 


S 3 ‘ d~ |a 2 (x ' 
4 


x P " “ ypi 
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The factor g is defined in (5.6), DFD in (5.11), and sin 6 in (5.8). The 
final expressions for the diffracted field in the umpriraed coordinate 
system are as follows: 

H d (P.) = xH d + yH d + zH d (5.31a) 

2 x J y z 


where 



if T lies on an elliptical cylinder; and 



) H d = -h„ cos 0, - h_ sin 0_ 
y 2 3 3 3 

H d = h- sin 9, - h, cos 0_ 
z 2 3 3 3 


(5.31b) 


(5.31c) 


if T lies on an elliptical cone. The corresponding diffracted electric 
field E d at ? 2 is calculated from 

E d (P 2 ) = 120n d~ 1 [H d (P 2 ) x d 4 ] (5.32) 

5.7. Detour Parameter 

The diffracted field H d calculated from (5.5) is not valid when obser- 
vation point P 2 is near the reflected shadow boundary which is defined by 
(5.25). To detect if P 2 is indeed so, we may calculate the so-called 
"detour parameter" of the reflected field (Section VI of [18]). 

5 = e[k(d 3 + d 4 - ,d l - d 2 )] 1/2 (5.33) 


Here e is the shadow indicator of the reflected field defined by 
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+1 , if is in the shadow side of the reflected field 


-1 , if P„ is in the lit side of the reflected field 
2 


It may be shown that 


e 


-c 


Sgn -cos 


2 


<f> + Q 1 ) 


(5.34) 


The square root in (5.33) takes positive real, negative imaginary, or 
zero value. When the caustic of the reflected field falls on the reflected 
shadow boundary, £ is imaginary; otherwise, C is always real. Detour 
parameter F, in (5.33) becomes zero when is exactly on the reflected 
shadow boundary, because the diffraction point 0 1 * and reflection point 
0 r coincide (Figure 2). Following the numerical study in [24], we take 

Ul - 2 (5.35) 


as on the on-set point. Thus, if )c| > 2, is considered to be away from 
the reflected shadow boundary, and the diffracted field if* in (5.5) is 
valid. If |c| < 2, P^ is considered to be near the reflected shadow boun- 
dary, and we must replace if* by if 0 as stated in (2.9). 

When F, is small, (d^ + d^) in (5.33) is nearly equal to (d^ + . In 

many practical problems, the reflection point 0 r and diffraction point 
0^ may not be determined with great precision. Thus, when F, is small, a 
direct computation of F, from (5.33) can have a numerical accuracy problem. 

To circumvent this possible problem, we have given below an alternative formula 
for F, when its value is small: 
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5 


cos 


( 4 * + 


4 1 ) 


_ /d. [1 + (d./R )] 

sin 0 /2k ■ -- 

/l + (d 2 /R[)/l + (d 2 /R2) 



(4> + 4 i ) 


__ (DF)/d 4 
sin 0 /2k ' ( DE )F) » 


if |M > 0 


(5.36) 


where (DF) of the reflected field is defined in (4.8), and (DDF) of the 
diffracted field in (5.11). The derivation of (5.36) is given in 
Appendix A of [14]. In all of the following computations, we use (5.33) if 
|5 | > 2, and (5.36) if |5 | < 2. 


5.8. Uniform Asymptotic Theory 

We shall calculate H° by the UAT developed in [14], [19], [21]. The for 
mula reads 


H D (P 2 ) = H d ( P 2 ) + 


F«) - F(5) " j (1 - e > 


H r (P 2 ) 


(5.37) 


where F is a Fresnel integral defined by 


\ “1/2 jir/4 r -jt^ . 

F(z ) = tr e J J e J dt 


(5.38a) 


and 


F(z) = — exp 

2z/tt 


• / 2 TT, 

"J (z + 4) 


(5.38b) 


The factor (1 - e)/2 in (5.37) is one if P 2 is in the lit region of the 
reflected field, and zero, if ? 2 is in the shadow. As expected, reduces 

■f(J 

to H when P 2 is away from the reflected shadow boundary. This is because 
of the fact that for |? | * =», 
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(5.39) 


[ F( O - F(C)1 ~j (1 - e) 

When is near or on the reflected shadow boundary, fT* in (5.37) is always 
finite, and compensates exactly for the discontinuity in ft® so that the 
total field it 1 " in (2.9) is everywhere continuous. 
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6. SECONDARY PATTERN COMPUTATION 


6.1. Aperture Field Theory 

Refer to Figure 3. From field equivalence principles, solutions for 
the far field may be obtained knowing the tangential fields at E . The 
aperture plane E is taken to be perpendicular to the z-axis. 

3 

Let us denote the tangential electric and magnetic fields at E . by 

3 

E and H , respectively. The field may be determined by using vector 
3 3 

potentials [25]. However, it is more convenient to express the far field 
directly in terms of the aperture fields. Let us define the following vec- 
tor quantities: 

l(u,v) = // E (x,y) e^ k ^ U5Cfvy ^dxdy (6.1) 

E 3 
a 

g(u,v) = // H (x,y) e'^ UX+Vy ^dxdy (6.2) 

E a 
a 

where 

u = sin 9 cos $ 
v = sin 9 sin <j> 
k = 2ir/X 

9 ,4> = spherical coordinates of far field point 

Since the aperture fields are tangent to z, let 
+ 

f = f x + f y 
x y J 

> * 
g = S x x + g y y 


(6.3) 

(6.4) 


(6.5) 

( 6 . 6 ) 
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From each of the 3 equivalence principle formulations, three different 
expressions arise for Eg and E^ of the far field, namely, 


(1) using E and H 

3 ct 


_ 1 ke 


-jkr 


0 4irr 


|^f x cos(j) + fySimji + Z q cos8 (g^cosiji - g^sirajoJ 


(6.7a) 


E 


-jkr 


J <(> 4irr 


J^cosB (f^cosiji - f^sirajO - Z o (gySin<j> + g^cos^ )j 

(6.7b) 


(2) using H 


jk Z e ~^ kr cos8 

E 0 = ^ (g y C0S * “ g x Sin<f,) 


(6.8a) 


jk Z q e Jkr 

% = (g y sln * + g x C0S4>) 


(6.8b) 


(3) using E 


„ ~j kr 


E = (f cosi}> + f sin<J>) 

9 2tt r x T y 


(6.9a) 


E = ^ ke _ — ° - S — (f cos$ - f sin$ ) 

$ 2nr y r x 


(6.9b) 


This method suits large apertures (in terms. of \) because the Fourier 

-V 

transforms of the aperture field, f and g, are highly peaked in the fre- 
quency domain. 
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This method is exact if the fields at were known everywhere. 
However, one must truncate E to finite dimensions in order to employ the 

cl 

FFT. 


6.2 Fast Fourier Transform 

To employ the FFT subroutine, the integrals (6.1) and (6.2) must be 
rearranged so that the form of the integral matches the definition given by 
the documentation of the FFT subroutine. For this particular 2-dimensional 
FFT subroutine, the function being considered is assumed to be periodic in 
x and y with period 1 in x and y. Hence, the aperture grid as shown in 
Figure 12 must be scaled accordingly. 

Many manipulations are required in order to use the FFT. These mani- 
pulations are carried out for f x . Similarly, this can be applied to obtain 
g x » and gy. From (6.1) and (6.5) 


f x (u,v) = 


y 2 x 2 

/ / E ax (x,y) 

y l X 1 


jk(ux+vy) 


dxdy 


( 6 . 10 ) 


where E is the x-component of £ (x,y) in (6.1). Using the substitutions 

3X 3. 


X_ = (x-x. )/(x 2 ~x. ) 


1= (y-y 1 )/(y 2 -y 1 ) ! 


( 6 . 11 ) 

( 6 . 12 ) 


(6.11) becomes 


1 1 


f x (u,v) = Kj / j E^ x (xj+ _X_(x 2 _ x j ) , y^ ^SV 2 ~ Fj) ) 


j'^( u ( x 2~ x l )— + v(y 2 -y l ) — ^ ) 

2 dXdY 


(6.13) 


where 
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K ! = (x 2 -x 1 )(y 2 -y 1 ) e 


jf^uXi+vyj) 


(6.14) 


The expression ^ ax ( x ^ + ( x 2 ~ x l © under the integral may be 
interpreted as E (x,y) scaled to F(a,b) where CKa<l and 0<b<l (see Figure 
12). F(a,b) may be approximated by 


« 2 m 2 


F(a.b) - I l C e J2*W + "I> , 
L mn ’ 


(6.15) 


n=N ^ m=M^ 


where C are the Fourier coefficients obtained using the FFT subroutine, 
mn 

For a 32x32 FFT, 


N x - Mj - -15 


N 2 = M 2 = 16 


From (6.13) and (6.15), 


1 1 N 2 M 2 (mX+u(x2“Xi))x + (nX+v(y2~yi))Xl 


, 2tt 


f (u,v) = K J / / X l c 

*»■ 1 A n XT \M 


0 0 n=N^ 


mn 


dXdY 
(6.16) 


Interchanging the summation and integral signs and noting that 


1 

/ e 
0 


sin(£f) 
dx - e 

TT C 

X 


TTC 

X 


(6.17) 


the expression for f^ is 


N 2 M 2 

f K (u,v ) - K J I C nn o 
n=N^ m=M^ 


j^{ntf-n) sin[^{mX+u(x 2 -xi)) ]sin[^{nX+v(y 2 -yi)) ] 
it (mX+u(x2 - xi )) x ( nX+v(y2~yi ) ) 


(6.18) 


47 



where 


Jf (u(x 1+x2 ) + v(yl+y2)) 

K - (x 2 -x 1 )(y 2 -y 1 )e (6.19) 

In summary, to evaluate the integral (6.10), we first calculate {C 1 of 

mn J 

(6.15) by FFT and then calculate f x (u,v) via (6.18). 

6.3 Gain Normalization 

The secondary pattern is usually decomposed into two orthogonal 
polarizations. Following Ludwig's definition 3 [26], the following unitary 


polarization vectors are introduced. 

R = 0(ae^cos<}> + bsin$ ) + $ (-ae-^sin<j> + bcosifi ) (6.20) 

C = 0(ae -^sin<}> - bcosij) ) + <}>(ae -^costj> + bsin<j>) (6.21) 

Let the secondary pattern be expressed as 
“jkr p. /> — » 

£ - S—. [_9 U 0 + $ U(j) J (6.22) 

The reference-polarization and cross-polarization expressions of $ are 

Reference-pol of E = (e»(R*)*) (6.23a) 

Cross-pol of E = (E«(C*)*) (6.23b) 


The second conjugate operation in (6.23) results from the change in direc- 
tion of the field after being reflected by the reflector. 

The directivity for the reference polarization is defined by 



Similarly, the directivity for the cross polarization is defined by 
4n |e* c | 2 /z 

D c (Q,4>) — f (6.25) 

inc 

where p ^ nc is the incident power radiated from the feed. Noting that 
e ^ r /r factor is common in both the secondary pattern (6.22) and the inci- 
dent radiated field (3.6) the directivity formulas are 


4u |u (ae^cosij) + bsirnji) + U (-ae^sin<j> + bcos4>)| 2 
D R (0,<j>) = - — p * (6.26) 


Ait | U (ae ~^sin<j> - bcos<}> ) + U (ae ^cosiji + bsin<j>)| 2 

D c O ,<{,) = - g — p < 6 ‘ 27 > 

o inc 

Thus, for any feed polarization (a,b,t), (6.26) and (6.27) give the 
reference and cross-polarization directivities. 
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7 . NUMERICAL EXAMPLES 


The present aperture integration method for calculating the secondary 
pattern of a reflector entails the following steps: 

(i) Input the reflector geometry, eq. (2.1), and the boundary 
description. 

(ii) Input the feed pattern as described in eq. (3.4). 

(iii) Calculate GO contribution to the aperture field. 

(iv) Calculate edge-diffraction contribution to the aperture field. 

Use either UAT [14] or UTD [15] for aperture grid points near the 
shadow boundary. 

(v) Use the FFT, eq. (6.15), to obtain the far field. 

(vi) Decompose the far field pattern into reference-pol and cross-pol 
components using eq. (6.23). 

(vii) Use eqs. (6.26) and (6.27) to obtain the reference and cross- 
polarization directivities. 

In this section, we shall present some numerical results to establish 
the numerical accuracy of the present method. Near field, far field and 
scan data for a large reflector are presented. 

7.1. Effects of Aperture Grid Size and Location on Secondary Pattern 

The aperture field theory used to determine the secondary pattern is exact 
if the tangential fields are known everywhere on the aperture plane £ . 

cl 

When employing the FFT, £ is truncated. To minimize the amount of com- 

cl 

puter time spent, the size of £ should be as small as possible while cap- 
turing almost all of the field. To this end, a study of varying the aper- 
ture grid size and location of £ fl was performed to determine their effects on 
the secondary pattern. 
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The reflector used for this study is an offset parabolic reflector 


(see Figure 13) with typical f/D p and f/D values: 

■§-- 0.40, 1.00 (7.1) 

P 

The corresponding values of 0^ , 0^ » ami H are 

Oj - 14.25° (7.2a) 

0 2 - 64.01° (7.2b) 

H - 0.25D (7.2c) 


The reflector is being illuminated by a y^-polarized feed at focus with 
10 dB feed taper. Assuming (cos 9)^ type patterns, eq. (3.5), the E- and 
H-plane feed patterns are given by 

U E (0) - (cos 9) 11,82 (7.3a) 

U H (0) - (cos 9) 11,82 (7.3b) 

Refer to Figure 14. Let the aperture plane Z be located a distance L 

3L 

away from the focal plane. The focal plane is located at z * f. Plots of 
the tangential components of the aperture field along two cuts of Z are 

cl 

obtained as the distance L is varied. These two cuts are 

x-cut: y=H+0.5D, z«f+L 

y-cut: x =* 0 , z - f + L 

The diameter of this test reflector is 50 X. Due to the choice of feed 
polarization, only the y-component of the electric field, E^, is plotted. 
Figures 15(a)-(c) are plots of |E | for a x-cut for L * 0, 10X, and 20X, 


51 




(b) x-cut 


Figure 14. Two cuts of aperture plane E a at distance L away from focal plane. 
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E-FIELD (V/m) 



x/X 


Figure 15b. | Ey j for a x-cut, L * 10 X 
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-FIELD (V/i 








respectively. Figures 16(a)-(c) are plots of |E^| for a y-cut for L - 0, 
10X, and 20X, respectively. 

From Figures 15(a)-(c) and 16(a)-(c), note that the aperture field 

does not spread as L increases. This is so because the feed is at the 

focus. To demonstrate that the aperture field does spread as L increases, 

let's move the feed toward the reflector along z f by an amount d f (see 

Figure 17). Figures 18(a)-(c) are plots of (E^l for a x-cut for d^ » 2X 

at L = 0, 10X, and 20X, respectively. Figures 19(a)-(c) are for the y-cut. 

Figures 18(a)-(c) and 19(a)-(c) show a slight increase of |E^| for grid 

points near the edge of the aperture grid. Also note that the width of the 

region where the GO field is nonzero increased when L is increased. Thus, 

to minimize the field strength for a fixed size of E fl , L = 0 was chosen. 

Let the aperture grid be W by W. Choosing the criterion that the field at 

the edge of E is at least 20 dB below the maximum field value on E , W = 

3 Si 

1.1 4D centered at the midpoint of the projected aperture was chosen. 

7.2 DBS Antenna 

A direct broadcast satellite (DBS) antenna was designed by Lee et al. 
[27]. Pattern computation programs using the Jacobi-Bessel series tech- 
nique [4] were developed by Y. Rahmat-Samii. A parabolic reflector was 
used and it is described by 

D =* dish diameter = 108.148X (7.4a) 

f = focal length =» 94.867X (7.4b) 

H = offset height = 16.865X (7.4c) 

The feed is located at the focus and the primary pattern, eq. (34), is 
described by 
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E-FIELD (V/m) 


0.020 



y/X 

Figure 16a. | Ey j for a y-cut, L « 0 
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E-FIELD (V/m) 



y/X 


Figure 16b. E y for a y-cut, L ■ 10 A 
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E-FIELD (V/m) 



Figure 16c. | E y | for a y-cut, L » 20X 
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E-FIELD (V/m) 



x /A 


Figure 18a. |Ey,| for a x-cut, =* 2X, L ■ 0. 
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E-FIELD (V/m) 



x /X 

Figure 18b. |E^| for a x-cut, * 2X, L ■ 10 X. 
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E-FIELD (V/m) 



Figure 18c. |Eyj for a x-cut, * 2 A, L - 20 A. 
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E-FIELD (V/m) 



y/X 

Figure 19a. |E y [ for a y-cut, * 2X, L ■ 0. 
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E-FIELD (V/ 







(7.5a) 


a = 1//2 , b » 1//2, ^ - 90° (RHCP) 

O E (0) = (cos 0) 3 * 6 (7.5b) 

U R (0) - (cos 0) 2 ' 8 (7.5c) 

Using this aperture integration technique, the secondary pattern for <j> » 0° 
was computed using GO and GTD constructions. As shown in Figure 20, two 
patterns are superimposed with Rahmat-Samii* s results. One pattern used GO 
fields only and in computing the other pattern, the edge-diffracted field 
was included in the aperture field calculation. The gain and sidelobe 
levels of these two patterns are tabulated versus Rahmat-Samii' s results in 
Table 3 below. 

TABLE 3 

COMPARISON WITH RAHMAT-SAMII' S RESULTS 



Rahmat-Samii 

GO + GTD 

GO 

main beam 

48.28 

48.33 

48.32 

1st sidelobe 

28.42 

28.29 

28.42 

2nd sidelobe 

22.29 

22.18 

22.93 

3rd sidelobe 

18.05 

18.02 

18.12 

4th sidelobe 

14.95 

14.96 

13.40 

5th sidelobe 

12.39 

11.85 

11.14 

6th sidelobe 

10.31 

9.05 

8.41 


Due to the limited amount of computer working space, a 48 x 48 FFT was 
used. Despite the coarse sampling, the results in Table 3 are in good 
agreement. 

The scan performance of this reflector has been studied by Hung [28]. 
The Fourier-Bessel series technique is used to compute the secondary pat- 
tern. Using the same feed, eq. (7.5), the feed was displaced 5.8 in the 
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REF-POL DIRECTIVITY (dB) 



9 (degrees) 


Figure 20 • Secondary pattern of a DBS antenna using Jacobi-Bessel series 
GO 9 and GTD • 



-x^. direction to produce a pattern with peak directivity at 0 = 3° (see 
Figure 21). The H-plane scan pattern was also computed using the aperture 
integration technique with and without the edge-diffracted field. Figure 
22 depicts the H-plane scan pattern using the aperture integration and 
Fourier-Bessel series techniques. The patterns are in good agreement. 


7.3. TRW Antenna 

A dual reflector antenna system was designed by TRW for NASA-Lewis 
Research Center. The dual reflector is an offset Cassegrain reflector. 


The main reflector is parabolic and it 

D = dish diameter = 257. 89X 
f = focal length = 318. 74X 
H = offset length ■ 135. 51X 

The reflector is being illuminated by 

18 dB edge taper. 

Two cuts of E at the focal plane 
a r 

x-cut: y = H + 0.5D, z = f 
y-cut: x = 0 , z = f. 


is described by 

(7.6a) 

(7.6b) 

(7.6c) 

. y^-polari^ed feed at focus with 
were taken. 


The magnitude and phase of the y-coraponent of the electric field, E 
plotted in Figures 23(a)-(b) for the x-cut and in Figures 24(a)-(b) 


, are 

y 

for the 


y-cut. 


The secondary pattern for <J> = 0° of this main reflector is shown in 
Figure 25. The key features of the ref erence-pol directivity plot are 
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REF-POL. DIRECTIVITY (dB) 





E-FIELD (V/m) 



x/X 


Figure 23a. 


Amplitude of E for a x-cut of TRW main reflector. 
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PHASE OF E y (deg) 


360 



x /X 


Figure 23b. 


Phase of E for a x-cut of TRW main reflector. 

y 
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20 . 160 . 200 . 240 . 280 . 320 , 

y/A 

Figure 24a. Amplitude of E for a y-cut of TRW 


, 360 . 400 . 

main reflector 





PHASE OF E y (de< 



Figure 24b. Phase of E for a y-cut of TRW main reflector 




REF-POL. DIRECTIVITY (dB) 



9 (degrees) 


Figure 25. Reference-pol. directivity plot of TRW main reflector, 4> = 0°. 
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Figure 26 is a plot of the cross-pol directivity. 

The hyperboloid subreflector has a magnification factor of 2. In 
approximating the performance of this dual reflector system, we employ the 
equivalent paraboloid [29]. The equivalent paraboloid is described by 

D = 257.89 (7.7a) 

r- 2(318.74) - 637.48 (7.7b) 

H - 135.51 (7.7c) 

The corresponding secondary pattern for =* 0° is shown in Figure 27. 
The gain, HPBW, and sidelobe level (SLL) are 



Figure 28 is a plot of the cross-pol directivity. 

Scan performances of the TRW main and equivalent reflector were 
studied. The scan plane chosen is the x^y^ plane and E- and H-plane scans 
were performed (see Figure 21). 

For E-plane scan, the coordinates of the feed are 
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CROSS-POL. DIRECTIVITY (dB) 


60 . 



Figure 26. Cross-pol. directivity plot of TRW main reflector, <f> ■ 0°. 
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REF-POL. DIRECTIVITY (dB) 



9 (degrees) 


Figure 27. Reference-pol. directivity plot of TRW equivalent reflector, <j> 




CROSS-POL. DIRECTIVITY (dB) 



6 (degrees) 


Figure 28. Cross-pol. directivity plot of TRW equivalent reflector, <j> ™ 0 




X = 0 


(7.8a) 


y * — f tan (N ♦HPBW) cos 9 £ 
z - f - f tan (N *HPBW) sin 9 c 

where 


(7.8b) 

(7.8c) 


f = focal length 

„ . -1 ,H + 0.5D X 

9 c " 2 tan ( 2f } 

N = number of beamwidth's scan 
HPBW =* beamwidth of secondary pattern ( <j> 
For H-plane scan, the feed coordinates are 
x ■ — f tan (N •HPBW) 
y =■ 0 
z ■ f 


- 90°) with 
given by 


(7.9) 


feed at focus. 


(7.10a) 

(7.10b) 

(7.10c) 


where 

f =* focal length 
N number of beamwidth' s scan 

HPBW » beamwidth of secondary pattern ( <j> * 0°) with feed at focus. 

As the feed moves away from the focus, the secondary pattern degrades. 
The degradation may be characterized by peak gain loss and half-power beam- 
width. These 2 f igures-of-merit are plotted versus number of beamwidths 
scanned for the following four cases: 

(1) TRW main reflector, E-plane scan (see Figures 29-30) 

(2) TRW equivalent reflector. E-plane scan (see Figures 31-32) 

(3) TRW main reflector, H-plane scan (see Figures 33-34) 

(4) TRW equivalent reflector, H-plane scan (see Figures 35-36). 

As expected, the scan performance of the equivalent reflector is much 
better than for the main reflector due to its larger f/D value. For the 
same number of beamwidths scanned, the peak gain loss and amount of beam 
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PEAK GAIN LOSS (dB) 



BEAMWIDTHS SCANNED 


Figure 29. Peak gain loss versus beamwidths scanned for E-plane scan of TRW 
main reflector. 




HALF-POWER BEAMWIDTH (deg) 



BEAMWIDTHS SCANNED 


Figure 30. Half-power beamwidth versus beamwidths scanned for E-plane scan 
of TRW main reflector. 




PEAK GAIN LOSS (dB) 



Figure 31. Peak gain loss versus beamwidths scanned for E-plane scan of TRW 
equivalent reflector. 
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HALF-POWER BEAMWIDTH (deg) 



Figure 32. Half-power beamwidch versus beamwidths scanned for E-plane scan 
of TRW equivalent reflector. 




PEAK GAIN LOSS (dB) 



Figure 33. Peak gain loss versus beamwidths scanned for H-plane scan of TRW 
main reflector. 




HALF-POWER BEAMWIDTH (deg) 



BEAMWIDTHS SCANNED 


Figure 34. Half-power beamwidth versus beamwidths scanned for H-plane scan 
of TRW main reflector. 
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PEAK GAIN LOSS (dB) 



Figure 35. Peak gain loss versus beamwidths scanned for H-plane scan of TRW 
equivalent reflector. 
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HALF-POWER BEAMWIDTH (deg) 



Figure 36. Half-power beamwidth versus beamwidths scanned for H-plane scan 
of TRW equivalent reflector. 




broadening were much less for the equivalent reflector for both E- and H- 
plane scans. Also, the cross-pol directivity was much lower for the equiva- 
lent reflector. Figures 37(a)-(b) and 38(a)-(b) are the ref erence-pol and 
cross-pol directivity plots for 6 beamwidth scan in the H-plane for the 
TRW main and equivalent reflectors, respectively. With respect to the peak 
directivity, the maximum cross-pol directivity value is -39.15 dB for the 
equivalent reflector and only -27.95 dB for the main reflector. 
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REF-POL. DIRECTIVITY (dB) 



9 (degrees) 


Figure 37a. Ref erence-pol. directivity plot of TRW main reflector for 6 
beamwidth scan in the H-plane. 
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9 (degr 


Figure 37b. Cross-pol. directivity plot 
width scan in the H-plane. 



TRW main reflector for 6 beam- 






A1IAI103; 





CROSS-POL. DIRECTIVITY (dB) 



6 (degrees) 

Figure 38b. Cross-pol. directivity plot of TRW equivalent reflector for 6 
beamwidth scan in the H-plane. 
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8 . CONCLUDING REMARKS 

We have developed a method of calculating the secondary pattern of an 
arbitrarily shaped reflector illuminated by a feed with arbitrary polariza- 
tion. An edge-diffracted field was added to the geometrical optics field in 
the aperture field calculation. By employing the FFT, the secondary pat- 
tern is computed very efficiently. The results for the secondary pattern 
are in good agreement with those obtained by the physical optics integral. 
Furthermore, this method can be conveniently extended to secondary pattern 
computation of multiple reflector systems, which will be done in the next 
phase of this project. 
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APPENDIX A 


INSTRUCTIONS TO RUN COMPUTER PROGRAMS 
To run the computer programs on the University of Illinois CYBER 175 
system to compute the secondary pattern of an offset parabolic reflector, the 
user should perform the following 

(1) create input file TAPE2 

> 

(2) run program BALI 

(3) run programs BPT, BFFT 

(4) for plot, use BPL0T3 

The computer commands corresponding to steps ( 1 )— (4) are 

(first create TAPE2) 

R 

BALI 
R 

P.LOAD(BPT,BFFT); EXECUTE 
R 

BPL0T3 

The secondary pattern is stored in TAPE13. 

The input file TAPE2 consists of 14 lines. These inputs are described 
below. 

1. H,D,f ,z Q — offset parabolic reflector parameters as depicted in 
Figure 13 

2. 1 — cylinder boundary or 2 — cone boundary 

3. x c , y c> , Kj- — cylinder boundary parameters or P, 0 i» 0 2» 0 3 — 
cone boundary parameters as shown in Figure 5 

4. Y 1 ,Y 2 > Y 3 — Eulerian angles (see Figure 7) 
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5 


* X 1 » y l » Z 1 — source point coordinates in main reflector coordinate 
system. 

6. a,b,i|; — polarization of feed 

7. q E ,q H — feed taper parameter for E- and H-plane primary patterns 

8. L — location of aperture plane beyond the focal plane as shown 
in Figure 14a (taken to be zero) 

9. type in 1 or 2 

x 1 ,x 2 .y 1 ,y 2 are bounds of E a (see Figure 12a), N is the 
number of points in sampling. 

12. 0 — use UAT or 1 — use UTD, 0 — compute GO and edge-diffracted 
field or 1 — compute GO field only 

13. 0p0 2 ,NFFP — For a theta cut, 0^ is the first value and 0 2 is the 

last value of theta. NFFP is the number of 0 values for the 

secondary pattern 

14. 1 — use E-fields only for the aperture field calculation or 2 — 
use H-fields only, PHI — this is the constant <|> value when 
computing the secondary pattern. 

For an arbitrarily shaped reflector, the user should use the binary ver- 
sion of MAIN instead of BALI. However, lines 97 to 107 must be replaced with 

either the coordinates of the points that define the reflector or the new 
equations that describes the reflector. For an arbitrary reflector, line 1 of 
input file TAPE2 has no meaning. Thus, delete line 41 of MAIN and the Input 
file is lines 2 through 14 of TAPE2 described above. 


1U. x j , x 2 , N 

11* y, .yo.N 
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APPENDIX B 


VIEWGRAPHS FOR A PAPER ENTITLED "CALCULATION OF NEAR-FIELD OF A REFLECTOR BY 
GTD" BY P. T. LAM, S. W.. LEE, AND R. ACOSTA PRESENTED IN INTERNATIONAL 
IEEE/AP-S SYMPOSIUM NATIONAL RADIO SCIENCE MEETING, JUNE 1984 

GTD Analysis of Re-f lector Antennas 


T.T. Law and S. U. Lee. 

Electromagnetics La!o 
University of Illinois 

and 

R • A CO$ta 

N as a.- Lewis feearch Center 


Outline. 

1. Troblem description 

2. Motivation 

3. Approach 

4. Results 

5 . Conclusion 
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Methods 
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Advantages or aperture integration (Al) 
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■■Held measurements 


3 . Accuracy is Comparable to pkqsfCftJ 
optics CTO). The use of -the, FFT 
makes it numerically eff-icieat- 
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What's new m +his AX approach? 


1. Ed^je- diffracted -field is included in 
'fhe. aperture Ae(d calculation. 

2. Surface of reflector is Complefeloj 
Arbitrary. 

3. Eic^e of reflector not restricted 

Ao a Circular cur ye, It may be an 
arbt'traryj Curye lyiVuj on an elliptical 
Co^e or cylinder. 
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CONCLUSION! 

& UJi'-tK dikes AI approach, -the- mam beam and 
Side. lobe returns of -the second anj pattern, 
i S j accurately computed. 

The ref lector shape is completely arbitrary. 

This -technique is numerically e-ff&ent ■ 
by employing the F FT. 

N^ar -field measurements may be checked, 

-)£- This AI -forw^lubon iS Suitable, -for 
extension ~to multiple reflectors . 
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-%■ extend technique -to 
dual ret lector System 

#■ array -feed 

-)f multiple reflectors 
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